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High Curie temperature of 900 K has been reported in Cr-doped AIN diluted magnetic semicon- 
ductors prepared by various methods, which is exciting for spintronic applications. It is believed that 
N defects play important roles in achieving the high temperature ferromagnetism in good samples. 
Motivated by these experimental advances, we use a full-potential density-functional-theory method 
and supercell approach to investigate N defects and their effects on ferromagnetism of (Al,Cr)N 
with N vacancies (Vk)- We investigate the structural and electronic properties of the Vn, single Cr 
atom, Cr-Cr atom pairs, Ct-Vn pairs, and so on. In each case, the most stable structure is obtained 
by comparing different atomic configurations optimized in terms of the total energy and the force 
on every atom, and then is used to calculate the defect formation energy and study the electronic 
structures. Our total energy calculations show that the nearest substitutional Cr-Cr pair with the 
two spins in parallel is the most favorable and the nearest Cr-ViM pair makes a stable complex. Our 
formation energies indicate that Vn regions can be formed spontaneously under N-poor condition 
because the minimal formation energy equals -0.23 eV or Cr-doped regions with high enough 
concentrations can be formed under N-rich condition because the Cr formation energy equals 0.04 
eV, and hence real Cr-doped AIN samples are formed by forming some Cr-doped regions and sepa- 
rated Vn regions and through subsequent atomic relaxation during annealing. Both of the single Cr 
atom and the N vacancy create filled electronic states in the semiconductor gap of AIN. N- vacancies 
enhance the ferromagnetism by adding to the Cr moment each, but reduce the ferromagnetic 
exchange constants between the spins in the nearest Cr-Cr pairs. These calculated results are in 
agreement with experimental observations and facts of real Cr-doped AIN samples and their syn- 
thesis. Our first-principles results are useful to elucidating the mechanism for the ferromagnetism 
and exploring high-performance Cr-doped AIN diluted magnetic semiconductors. 



I. INTRODUCTION 

The diluted magnetic semiconductor (DMS) is of great 
interest because it has merits of both semiconductor 
and spin magnetism. Especially the possibility of high- 
temperature ferromagnetism in DMSs has inspired sci- 
entists from all the world to explore new promising 
DMS materials and study their physical properties for 
spintronic applicationsi"^. (Ga,Mn)As has been stud- 
ied substantially because of its close relation with im- 
portant GaAs semiconductor, and its ferromagnetism 
is believed to be induced by hole carriers^. Gen- 
erally speaking, it is believed that at least 500 K as 
Curie temperature is necessary for real spintronic de- 
vices to work at room-temperature^. Curie tempera- 
ture higher than 500 K is achieved in various transition- 
metal doped oxides^ii^iSiiOiiiii^iHiii and nitridesi^ii^. 
For Cr doped AIN, or (Al,Cr)N, high Curie tempera- 
ture of 900 K has been reported by several groupai^'i^'i^. 
It is shown that the high temperature ferromagnetism 
in (Al,Cr)N materials is intrinsic, not produced by sec- 
ondary phase o^^'^'^i^^'^^i^°'^^i^^ . Naturally, such high 
Curie temperatures make the materials promising for 
practical spintronic applications^. Although ferromag- 
netism of (Ga,Mn)As DMS materials is believed to be 
induced by their effective carriers^ii., the high temper- 
ature ferromagnetism of (Al,Cr)N materials must be 
caused by different mechanism because some single-phase 
(Al,Cr)N samples with high Curie temperatures are 



highly insulating^i^ (the same as pure AIN-^) or in the 
regime of variable-range hopping conduction24 (similar to 
(Ga,Cr)N^-). The ferromagnetism should be attributed 
to pd hybridization of Cr d states and N p one o^^i^'''^^ . 
On the other hand, defects and impurities are believed to 
play important roles in DMS materials of wide semicon- 
ductor gap a-^'*'^^i^^'^^i'^°i'^^'^^ '^. Nitrogen pressure during 
synthesis is shown to affect the ferromagnetism in re- 
sultant sample o^^i^^'^" . It appears that reduced N con- 
tent enhances the ferromagnetism, being in contrast with 
first-principles-calculated trend of N- vacancies in the case 
of GaN-based DMSs^. Therefore, it is highly desirable 
to investigate effect of N content on the high temperature 
ferromagnetism in (Al,Cr)N materials, and in return its 
elucidation is helpful to obtain deeper insight into the 
essential mechanism of the high temperature ferromag- 
netism. 

In this paper we use a full-potential density-functional- 
theory method and 3x3x2 supercell approach to inves- 
tigate N defects and their effects on ferromagnetism of 
(Al,Cr)N with N vacancies (Vn), namely Cra;Ali_a;Ni_j,. 
We obtain the most stable structures of AI36N35 (one 
in the supercell), CrAlssNsg (one substitutional Cr in the 
supercell), CrAlssNss (one Vn and one substitutional Cr 
atom in the supercell), Cr2Al34N3g (two substitutional 
Cr atoms in the supercell) , and so on by comparing their 
different atomic configurations optimized in terms of the 
total energy and the force on every atom, and thereby 
calculate their defect formation energies. Then we in- 
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vestigate the structural and electronic properties of the 
V^, single Cr atom, Cr-Cr pairs, and Ci-Vys pairs. With- 
out N vacancy, the nearest substitutional Cr-Cr pair with 
their moments orienting in parallel is the most favorable 
in total energy. The nearest pair of an N vacancy and a 
single Cr atom forms a stable complex of Cr^ and V^^ . 
Our calculated formation energies indicate that Vn re- 
gions can be formed spontaneously under N-poor condi- 
tion or high enough Cr concentrations can be formed un- 
der N-rich condition, and hence real Ali_a;Cr2.Ni_y sam- 
ples are formed by forming Cr-doped regions and some 
separated Vn regions and through subsequent atomic re- 
laxation during annealing. The doped Cr atom and the 
N- vacancy create filled electronic states in the semicon- 
ductor gap of AIN. Vn enhances the ferromagnetism by 
giving one electron to a Cr ion to add lyus to the Cr mo- 
ment, which is in agreement with experimental observa- 
tion that low N pressure during synthesis is helpful to ob- 
taining better ferromagnetic Ali.^-Cr^jNi.j, samplesSSi^. 
On the other hand, reduces the ferromagnetic ex- 
change constants between the spins in the nearest Cr- 
Cr pairs, which is consistent with experimental fact that 
too many N vacancies are harmful to the ferromagnetism 
when the synthesis temperature is high enough to allow 
doped Cr atoms to diffuse. These first-principles results 
are useful to exploring high-performance (Al,Cr)N DMS 
samples and understanding the high temperature ferro- 
magnetism. 

The remaining part of this paper is organized as fol- 
lows. We shall present our computational method and 
parameters in next section. In the third section, we shall 
present our main results of optimized magnetic structures 
(through geometric and internal atomic position opti- 
mization) and defect formation energies. In the fourth 
section, we shall present main electronic structures, in- 
cluding energy bands and electron density distributions. 
In the fifth section, we shall discuss the ferromagnetism 
correlated with N defects in real samples. Finally, we 
shall give our main conclusion in the sixth section. 



II. COMPUTATIONAL METHOD AND 
PARAMETERS 

We use a first-principles method and supercell ap- 
proach to study wurtzite AIN with doped Cr atoms and 
N defects (including nitrogen vacancies and nitrogen in- 
terstitials) . We use a 3x3x2 supercell as showed in 
Fig. 1 and substitute a Cr atom for an Al atom to 
model Cr doping concentration 2.78% realized in MBE 
samplesi^ii2iii^'2i. This supercell is not very large to sim- 
ulate real experimental conditions, but it is large enough 
to model the representative concentrations of both doped 
Cr and and main properties of Cr-Cr and Cr-Vw in- 
teractions, especially considering that our calculations 
are done with full-potential linear augmented plane wave 
(FLAPW) method within the density-functional theory 
(DFT)^^. Using such supercell is also supported by the 




FIG. 1: (color online). Schematic structure of Cr-substituted 
3x3x2 wurtzite AIN supercell, CrAlssNae. The centering 
largest (cyan or white) ball is the first Cr atom, the small- 
est (yellow or light gray) ones are N atoms, and the other 
(red or gray) ones are Al atoms. An N vacancy can be cre- 
ated by removing an N atom, and the second substitutional 
Cr atom can be introduced by substituting a Cr for an Al 
atom. For examples, Vni (or Vi) is made by removing Nl, 
and Cr3 by substituting a Cr atom for A13. 



following calculated results that the total energy of both 
CrAi-CrAi and CrAi-VN pairs substantially increases with 
increasing the CrAi-CrAi and CrAi-V^ distances, respec- 
tively. The supercell consists of one Cr atom, 35 Al 
atoms, and 36 N atoms and its space group is No. 156. 
Such a supercell for the pure AIN consists of 36 Al and 36 
N atoms. The second Cr atom is introduced by substi- 
tuting for one of Al atoms in the supercell. For examples, 
Crl is made by substituting a Cr atom for the Al atom 
labelled as All in Fig. 1. A nitrogen vacancy Vn is cre- 
ated by removing one N atom. For examples, Vni (or 
VI) is made by removing the N atom labelled as Nl. An 
interstitial N atom is introduced by adding an N atom 
in the supercell. We obtain an octahedral N interstitial 
Nt if adding an N atom at an octahedral site, or obtain 
tetrahedral N interstitial Nt if adding an N atom at a 
tetrahedral site. We optimize lattice constants to min- 
imize the total energy, and at the same time optimize 
internal atomic coordinates to make the force on every 
atom less than 0.04 eV/A. All the electronic structures 
and magnetic moments are calculated with equilibrium 
lattice constants. 

All the presented results are calculated with the full- 
potential linear augmented plane wave method within the 
density-functional theoryS^, as implemented in the pack- 
age WIEN2k^, although we also use pseudo potential 
method to do some preliminary structural optimizations. 
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The generalized gradient approximation (GGA)^ is used 
as the exchange-correlation potential for all our presented 
results, and some important cases are confirmed by local- 
density-approximate calculations—. Full relativistic cal- 
culations are done for core states and the scalar approx- 
imation is used for the others, with the spin-orbit cou- 
pling being neglected^. The radii of muffin-tin spheres, 
i?nit, are adjusted so as to achieve the best accuracy 
of self-consistent calculations. We use the parameters 
i?nit-Ji^max=8.0 and make the sphere harmonic expansion 
upto Zmax=10 in the spheres. The Brillouin zone integra- 
tions are performed with the special k-point method over 
a 5x5x3 Monkhorst-Pack mesh. A convergence stan- 
dard, / \pn — Pn-i\dr < 0.00005 (p„ and pn-i are output 
and input charge density, respectively) , is used for all our 
self-consistent calculations. 
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III. OPTIMIZED STRUCTURES AND DEFECT 
FORMATION ENERGIES 

We optimize fully the lattice structures and internal 
atomic positions for N vacancy (AlagNag + VN^AlagNas), 
Cr doped AIN (CrAlssNsg), Cr-atom plus N- vacancy 
(CrAlssNae -I- VN=CrAl35N35), Cr-atom plus tetrahedral 
N interstitial (CrAlasNae -I- Nt), and Cr-atom plus oc- 
tahedral N interstitial (CrAlssNse -I- No). We can have 
various atomic configurations for pairs of Cr-^N, Cr-Nx, 
and Cr-No with given distance, and they correspond to 
different total energies. As for CrAlasNsg -I- Vn, the to- 
tal energy changes with the distance between Cr and V^, 
and the lowest total energy is reached when Cr and 
are the nearest neighbors Gx-V\ (connected with an ap- 
proximately horizontal Cr-N bond). The next lowest to- 
tal energy (0.23 eV) is reached by Gx-V2. For double 
Cr doped system, the total energy depends on the Cr-Cr 
distance, and the lowest total energy is reached when the 
two Cr (Cr and Crl) are in the same layer and connected 
with the same N (consistent with previous resull4^). The 
next lowest total energy (0.28 eV) is reached by Cr-Cr3 
and other equivalent Cr-Cr pairs. These tends are shown 
in Fig. 2. 

TABLE I: The calculated equilibrium lattice constants (a/c), 
the relative volume change (^i;) with respect to the pure AIN, 
the magnetic moment per magnetic atom (M), and the elec- 
tronic density of states at the Fermi level (-Dp). The pure 
AIN phase is presented for comparison. 



Name 


a/c (A) 


&v (%) 


M (^ib) 


Df 


AlseNse+l^j 


9.65/10.29 


1.74 


1.000 


>0 


CrAl35N36 


9.59/10.26 


0.21 


3.000 


>0 


CrAlasNse-hVNi 


9.58/10.22 


-0.16 


4.000 





CrAlasNae-fNTi 


9.66/10.26 


1.28 


2.000 





CrAl35N36+Noi 


9.61/10.29 


1.15 


4.000 





AI36N36 


9.58/10.27 


0.00 


0.000 






Our systematic calculations show that a doped Cr 



FIG. 2: Relative total energy of the Cr doped AIN with N va- 
cancies for different Cr-VN distances (filled squares) and that 
of Cr doped AIN for different Cr-Cr distance (filled circles). 
Here Cr-Crn means the pair of the first Cr and the second Cr 
at the Aln site, and Cr-Vn the pair of the first Cr and the 
N vacancy at the Nn site, as are labelled in Fig. 1. The two 
smallest total energies, corresponding to the nearest Cr-VN 
and Cr-Cr pairs, are set to zero. 



atom contributes 3fiB to the total magnetic moment, and 
a neighboring N- vacancy adds IpB because its electron 
just below the Fermi level is transferred to the Cr d state. 
Actually, an isolate N- vacancy, even without any Cr dop- 
ing, has a moment of l/is, but N vacancies are in para- 
magnetic order because inter-vacancy spin interaction is 
equivalent to zero within computational error. As for in- 
terstitial N atoms, a tetrahedral N carries a moment of 
-IpB and an octahedral N has l/is, but they are too high 
in total energy to exist in real samples. Two neighboring 
Cr atoms, such as Cr-Crl and Cr-Cr3, have 6fiB because 
their ferromagnetic couplings arc favorable in energy. In 
fact, the ferromagnetic orientation is lower in total energy 
by 464 meV and 195 meV than the corresponding anti- 
ferromagnetic one for Cr-Crl and Cr-Cr3, respectively. 
The energy difference 464 meV is decreased to 163 meV 
when one N vacancy is at the nearest site, and is fur- 
ther reduced to 25 meV when two N vacancies are in the 
nearest neighborhood. This trend is consistent with ef- 
fect of N vacancy on Mn-Mn spin exchange constants—. 
The lattice constants, the relative volume changes, the 
magnetic moments per supercell of the most stable con- 
figurations are summarized in Table 1. The most stable 
is the N interstitial Nti (Nqi), the nearest to the Cr 
atom, in the CrAlasNse-l-NTi (CrAlssNsg-l-Noi), among 
all the tetrahedral (octahedral) N interstitials. The cor- 
responding results for the pure AIN are also presented for 
comparison. The DOS at the Fermi level. Dp, is zero in 
all the cases except for the CrAlasNae whose Fermi level 
is at the Cr d impurity bands and the AlsgNag-l-VN whose 
Fermi level is at the N-vacancy-induced in-gap band. 
We further calculate formation energy for all the cases. 
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The formation energy of the supercell including nv^ N va- 
cancies and ncr substitutional Cr atoms is defined as^ldl 

where £^(?T.yj^, ncr) is the total energy of the supercell, 
and (/xn, /^Cr, Mai) are chemical potentials of N, Cr, and 
Al. Using nf^^, /^Ai"^' ^^'^ f^Cr^ chemical poten- 

tials of N gas, Cr bulk, and Al bulk, we have /in — 
mIT'' + Mn, A^Ai = Mai''' + Mai, and Mcr = Mcr"*' + Mcr- 
The thermodynamical stability of AIN requires 

Mn + MAI — main- (2) 

The fact that there is no spontaneous formation of CrN 
implies 

MN + MCr < MCrN- (3) 

We can express the minimal formation energy as 

-fnv,Mr-«a(M^r''-Mir) (4) 
-Knv^„M5?--nc.(M'cT-MS'") 

where the superscripts max (min) means the maximum 
(minimum) of the chemical potential. The chemical po- 
tentials Mn'''' Mai"'' ^^'^ Mcr"' '^^^ ^'^ calculated in terms 
of total energies of N gas, Al bulk, and Cr bulk. The max- 
imal and minimal chemical potentials can be obtained in 
terms of equality ([2]) and inequality ([3]) and, therefore, 
the minimal formation energy can be calculated in terms 
of Eq. as long as McrN and main are known. 

TABLE II: The formation energies (AE) of a substitutional 
Cr atom and an N vacancy, and their pairs in AIN. The '4-', 
'0', and '-' superscripts mean that the extra charge is equiva- 
lent to e, 0, and -e, where e is the absolute value of the elec- 
tronic charge. The last row corresponds to the case that the 
distance between Cr and I^n is the longest within the 3x3x2 
supercell, and the fourth row the case that the distance is the 
shortest. 



Name 


AE (eV) 


Al36N36-|-l^Nl 


-0.23 


Al36N36-fl^N^ 


0.1" 


Cr"Al35N36 


0.04 


Cr-Al35N36 + V^l|l 


0.77 


Cr-Al35N36-HVN'5 


2.47 



Our calculated result is main = —3-33 eV for AIN, and 
we use known data McrN = —1.30 eV for CrN-^. For a 
single neutral N-vacancy in AIN under poor N condition, 
we have ny„ — 1 and ricr — 0, and thus the formation 
energy is -0.23 eV, while the formation energy of a single 
N+ vacancy is 0.1 eV— . For a neutral Cr atom under rich 
N condition, we have nv^ = and ncr = 1, and thus the 
formation energy is 0.04 eV. This small formation energy 



0.04 eV is reasonable and consistent with experimental 
high Cr doping concentration&i^iili2iiS2ii2£. A previous 
unreasonably high formation energy 3.01 eV— is in con- 
tradiction with the experimental high Cr doping concen- 
trations. Charged Cr atoms have much higher formation 
energies. For the nearest Cr-Vwi pair, we have ny„ = 1 
and ncr = 1, and thus the formation energy is 0.77 eV, 
while more separated Cr-T^^s pairs have higher forma- 
tion energies. These results are summarized in Table 
II. Interstitial N atoms, both tetrahedral and octahedral, 
have very high total energies and formation energies, and 
therefore are not presented in Table II. This implies that 
either dilute N vacancies or dilute Cr atoms are easy to be 
realized in AIN because the interactions can be neglected 
in the dilute doping limit. 



IV. ELECTRONIC STRUCTURES 

We present spin-dependent DOS of the AI36N35 (in- 
cluding one N vacancy in the supercell), the CrAlasNse 
(including one Cr substitutional at an Al site), and the 
CrAlssNas (including one Cr-t-Vwi complex) in Fig. 3. 
The valence bands for every spin channel are filled with 
108 electrons for the AIN supercell, and in total we have 
216 electrons in the valence bands. The absence of one 
N atom removes 3 p electrons and re-organizes the band 
structures, as shown in Fig. 3(a). Now there are three 
electrons in the gap of AIN and the valence band mani- 
fold is filled with 210 electrons, 105 for each spin channel. 
The N vacancy cannot donate electron to the conduction 
bands because the Fermi level is too low with respect 
to the conduction bands. This is different from the N- 
vacancy of GaN^. Fig. 3.(b) shows the DOS of the 
CrAlssNag. In this case we have 108 electrons in the va- 
lence band manifold for each spin channel and the Cr d 
states are in the gap of AIN. The d states are split into 
a doublet (labelled as TI') and a triplet (labelled as T). 
The former is filled with two electrons and the latter by 
one. The Fermi level is in the bands of the triplet. Hence 
one Cr contributes 3ms to the moment. In Fig. 3(c) we 
show the density of states of the CrAlasNas. The inter- 
action of the doped Cr and the N vacancy makes the 
filled singlet originating from the N vacancy move lower 
in energy and become two spin-split parts: 'a' and 'b', 
and at the same time makes the Cr d triplet split into a 
doublet (also labelled as T in Fig. 3(c)) and a singlet. 
The valence band manifold below the part 'b' consists of 
105+105 filled states. As a result, one Cr-|-VN complex 
contributes 4/is to the moment. 

We present in Fig. 4 electron density distribution of 
the N-vacancy-induced filled states, the parts 'a' and 'b' 
indicated in Fig. 3(c), of the CrAlasNas. The three 
planes are defined by three-site series: Cr-14i-A13, Cr- 
Vn-A12, and A12-yN-A13. Here Cr, All, A12, A13 are the 
atoms labelled with the same signs in Fig. 1. It can 
be seen that the charge density distribution at the N va- 
cancy is about 0.014e/a.u."^, where e is the absolute value 
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FIG. 3: (color online). Spin-dependent total DOS (black solid 
lines) and Cr d partial DOS (red or gray dash lines) of the 
AI36N35 (including one Vn) (a), the CrAlssNae (including one 
Cr atom) (b), and the CrAlasNas (including one Cr+VNi com- 
plex) (c). The upper half for each panel is for majority spin 
channel and the other for minority spin. For the AI36N35, 
the filled part just below the Fermi level contains on electron, 
and thus one Vn contributes I/ib. For the CrAlssNse, there is 
one electron in the small part T and two in the filled part of 
'IF, and hence the total moment is 3^b. For the CrAlasNss, 
both T and 'IF are filled with two electrons, and the total 
moment is 4/i_B. The small parts 'a' and 'b', both filled with 
one electron, are due to Vn. 



of the electronic charge. The total charge within the N 
vacancy is estimated to be approximately equivalent to 
l.le. The remaining 0.9e charge of the parts 'a' and 'b' is 
distributed at the nearest Cr site (0.6e) and N site (0.3e). 
This indicates that the N-vacancy state is substantially 
hybridized with the Cr d states and substantial charge is 
transferred to the vacancy from the neighboring Al (or 
Cr) atoms although there is no N atom at the site. 

We present in Fig. 5 the Cr d electron density dis- 
tribution of the CrAlasNae and the CrAlasNas. The 
three planes are the same as those in Fig. 4. For the 
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FIG. 4: (color online). Electron density distributions of the 
'a' and 'b' bands of the CrAlasNss, defined in panel (c) of 
Fig. 3. The al, a2, and a3 (bl, b2, and b3) show the electron 
density distributions of the 'a' ('b') bands in the three planes 
defined by Cr-VNi-A13, Cr-VNi-A12, and A12-T/ni-A13, where 
Cr, All, A12, A13 means the atoms labelled with the same 
signs in Fig. 1 and the Vni is at the Nl site, labelled as 
'-f ', connecting the Cr atom and three Al ones. The outmost 
(green or light gray) line represents 0.001 e/a.u."^. The density 
contour increment is 0.001 e/a.u.''. 



CrAlssNse, it is estimated approximately that 80% of 
the 'ir charge is around the Cr atom and 70% of the 
T charge around the Cr, and the Cr electronic configu- 
ration is . For the CrAl35N35, approximately 90% of 
the TF charge is around the Cr atom and 85% of the 'F 
charge around the Cr. The remaining charge is around 
the nearest N sites. Therefore, we can conclude that Cr 
electronic configuration in a Cr+^N complex is c?^, not 
d^ . The extra electron comes from the N vacancy. Thus 
the valent state of the Cr-|-VN complex is Cr~-|-V[^. 



V. FERROMAGNETISM CORRELATED WITH 
N DEFECTS IN REAL SAMPLES 

Because the minimal formation energy for forming a 
neutral N vacancy in AIN is -0.23 eV, a small concentra- 
tion of N vacancies can be formed spontaneously in AIN 
under poor N condition. An N vacancy contributes l/is 
to the moment, but our calculations show that spin inter- 
action between N-vacancy spins is estimated to be zero 
within computational error. Hence N vacancies alone 
cannot form ferromagnetism. Because the minimal for- 
mation energy for forming a substitutional Cr in AIN 
is 0.04 eV, dilute Cr doping alone is easy to be real- 
ized in AIN at room temperature under N-rich condition. 
High concentrations of Cr atoms have been doped into 
^^]\ ji6, 17^24, 29^30 ^ although the Cr formation energy may 
be larger when Cr-Cr interaction cannot be neglected. 
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FIG. 5: (color online). Electron density distributions of the 
parts T (al, a2, aS) and '11' (bl, b2, b3) of the CrAlasNae, 
defined in panel (b) of Fig. 3, and those of the parts T (cl, 
c2, c3) and 'IF (dl, d2, d3) of the CrAlssNas, defined in 
panel (c) of Fig. 3, where the three planes are the same as 
those in Fig. 4. It should be pointed out that the part 'F 
of the CrAlasNss consists of two electrons, but that of the 
other consists of one electron only. The outmost (green or 
light gray) line represents 0.001 e/a.u.''. The density contour 
increment is 0.001 e/a.u."^. 



This implies that there is no barrier of formation en- 
ergy in the dilute doping hmit. The Cr doped AIN can 
be ferromagnetic through some mechanism based on hy- 
bridization between local Cr-d states and extended N-p 
states^^^i^s. 

The minimal formation energy for forming a Cr-|-VN 
complex is 0.77 eV. As a result, the concentration of such 
complex is at most 2.3xl0~® per Al site when the syn- 
thesis temperature is at 600 K, becomes 10"'^ at 1140 K, 
and reaches to 1% at 1700 K. Each of the Cr-|-VN com- 
plexes contributes 4/is to the total moment. It should 
be pointed out that the maximal equilibrium concentra- 



tion of such complexes at 1000 K is lower at least by 
two orders of magnitude than the Cr concentrations of 
MBE-grown samples (nominally 2 ~ and 
other ones (nominally 3 - 35%)i^iiii^i29i30 (Al,Cr)N 
diluted magnetic semiconductors. Thus the synthesis of 
real samples is not through growth of the Cr+VN com- 
plexes. 
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FIG. 6: Schematic DOS of the AlaeNse (a), the AI36N35 (in- 
cluding one Vn) (b), the CrAlasNss (including one Cr-t-Vki 
complex) (c), the CrAlssNse (including one substitutional Cr 
atom) (d), and the Cr2Al34N3g (including two nearest substi- 
tutional Cr atoms) (e). The two numbers for the shadow re- 
gion denote the numbers of the majority-spin p valence bands 
and minority-spin ones. A solid dot means a filled electron in 
the bands. 

It is reasonable to suppose that there coexist N vacan- 
cies, Cr atoms, and Cr-f-^N complexes in real samples and 
the fractional concentrations are dependent on synthesis 
conditionslSilS. Because of N pressure fluctuation dur- 
ing synthesis, a substitutional Cr and an N vacancy can 
be formed easily in different regions. This implies that 
some N-vacancy domains and some Cr-doping domains 
can be formed simultaneously in separated different re- 
gions. The former are paramagnetic and the latter are 
ferromagnetic, but these domains are subject to ther- 
modynamical equilibrium of electrons in real samples. 
We present the schematic DOS and the corresponding 
Fermi levels of the AlsgNag, the AI36N35 (including one 
N vacancy in the supercell), the CrAl35N35 (including 
one Cr-l-l^i complex), the CrAl35N36 (including one Cr 
substitutional at an Al site), and the Cr2Al34N36 (includ- 
ing Cr and Crl as defined in Fig. 1) in Fig. 6 (a ~ e). 
The Fermi level of the CrAl35N36 is consistent with that 
of the Cr2Al34N36. Because the Fermi level of the N va- 
cancy is substantially higher than that of the Cr atom. 
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the electron of the N vacancy must transfer to the Cr site 
and this electron transfer results in pairs of and Cr^. 
The Coulomb attraction between and Cr~ will make 
them tend to move together and form the neutral Cr+^N 
complex. Both of the Fermi levels of the Vj^ and the neu- 
tral Cr+^N complex are consistent with that of Cr atom. 
The electron of the N vacancy tends to be transferred to 
the Cr atom and the neutral Cr+VN complex. 

The N vacancy is in favor of the ferromagnetism 
through forming the Cr+Vw complex which contributes 
AfiB to the moment and Cr-|-2VN complex which con- 
tributes bfiB- On the other hand, the N vacancy can be 
harmful to the ferromagnetism through destroying the 
ferromagnetic spin coupling between two nearest substi- 
tutional Cr atoms. We can figure out what will happen 
in terms of the two characteristic synthesis temperatures, 
Tsi and Ts2 (> Tsi), for Cra;Ali_:cNi_j,, where x is Cr 
concentration and y is N- vacancy concentration. We sup- 
pose that there is no second phase and hence set y < 2x 
in the following2£. When T < Tgi, neither N- vacancy 
nor Cr atom can diffuse. A little of the N-vacancy elec- 
trons will transfer to Cr domains but further electron 
transferring will be blocked by the built-in electric field 
formed by the transferred electrons. As a result, the 
main feathers of the domain patterns will remain. The 
Cr moments are coupled by thermally-activated variable- 
range hopping of electrons^, and thus the N-vacancy en- 
hances the ferromagnetism22i^. When Tgi < T < Ts2, 
N-vacancy can diffuse, but Cr atom cannot. All the N- 
vacancies and their electrons will transfer to Cr domains, 
and many Cr-I-Vw and Cr-t-2VN complexes are formed. 
In this case, the coupling of moments is the same, but 
the conductivity becomes weaker because an additional 
barrier must be overcome for the electron to hop2^, and 
thus the enhancement decreases with increasing y. When 
T > Ts2, both N-vacancy and Cr atom can diffuse. Cr 
atoms tend to form clusters-^ (similar to Cr- and Mn- 
doped GaN"'^''*^''*^) and the N vacancies tend to change 
the inter-Cr ferromagnetic coupling into antiferromag- 
netic coupling. In this temperature range, the N-vacancy 
is harmful to the ferromagnetism. However, The N va- 
cancy can be made to enhance the ferromagnetism by 
keeping the synthesis temperature below Ts2 so that Cr 
clusters are unlikely. This is in good agreement with 
experimental observation that N vacancies enhance the 
ferromagnetism in samples made at low temperatur o^^'^° . 

VI. CONCLUSION 

We have investigated nitrogen defects and their effects 
on the ferromagnetism of Cr-doped AIN with N vacancies 
using the full-potential density- functional-theory method 



and supercell approach. We investigate the structural 
and electronic properties of the Vn, single Cr atom, Cr- 
Cr atom pairs, Cr-V^N pairs, and so on. In each case, the 
most stable structure is obtained by comparing differ- 
ent atomic configurations optimized in terms of the total 
energy and the force on every atom, and then is used 
to calculate the defect formation energy and study the 
electronic structures. Without N vacancy, two Cr atoms 
have the lowest energy when they stay at the nearest Al 
sites and their moments orient in parallel. An N vacancy 
and a Cr atom tend to form a stable complex of Cr~ and 
and the moment of 4/iB origins from the Cr elec- 
trons. Our calculated formation energies indicate that 
Vn regions can be formed spontaneously under N-poor 
condition or high enough concentrations of Cr atoms can 
be formed under N-rich condition, but only tiny concen- 
trations of Cr-|-T44 complexes can be realized under usual 
experimental conditions. Both of the N vacancy and Cr 
atom create states in the semiconductor gap of AIN and 
the highest filled state of the N vacancy is higher than the 
filled Cr states. The enhances the ferromagnetism by 
giving one electron to Cr ion to add l/is to the Cr mo- 
ment, but reduces the ferromagnetic exchange constants 
between the spins in the nearest Cr-Cr pairs. 

The formation energy results imply that real Cr- 
doped samples with N vacancies are formed by form- 
ing Cr-doped regions and some separated Vn regions 
and through subsequent atomic relaxation during an- 
nealing. The enhancement of the ferromagnetism is 
in good agreement with experimental observation that 
low nitrogen pressure enhances the ferromagnetism of 
Ali_a;Crj;Ni_j^ when the synthesis temperature is low so 
that Cr atoms cannot diffuse^^i^. On the other hand, 
too many N vacancies are harmful to the ferromagnetism 
when the synthesis temperature is high enough to make 
Cr atoms diffuse and tend to become clusters of sec- 
ondary phases in which Cr-Cr spin interactions mainly 
are antiferromagneti o^^'^^i^^'^°i'^^ . However, we can make 
N-vacancies enhance the ferromagnetism by keeping the 
synthesis temperature low enough to avoid Cr diffusion. 
All these first-principles results are useful to exploring 
high-performance Cr-doped AIN DMS materials and un- 
derstanding their high temperature ferromagnetism. 
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